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THE ROLE OF MULTIPLE SCATTERING IN 
ONE -DIMENSIONAL RADIATIVE TRANSFER 


David Adamson 
Langley Research Center 

SUMMARY 

The usual methods of solving the radiative transfer equation yield answers which 
embrace all orders of scattering and thus shed little light on the underlying physical 
process. The present analysis examines the contributions of the various orders of 
scattering to the one-dimensional transfer of radiation. In the one -dimensional case an 
exact analytical solution exists and the problem reduces to that of expanding this exact 
solution in powers of the albedo for single scattering. Formulas are given which permit 
the calculation of any order of scattering in an atmosphere of arbitrary optical thickness, 
particle albedo, and asymmetry parameter. The results should aid in identifying those 
physical situations where only the lowest orders of scattering play a significant role and 
where appropriate approximate methods might provide results of acceptable accuracy. 

INTRODUCTION 

If one seeks to determine the radiation reflected from or transmitted through a 
plane parallel homogeneous atmosphere (slab) by solving, in the usual way, the radiative 
transfer equation, the answer one obtains embraces all orders of scattering. Ordinarily, 
this is all that is needed; however, such an approach provides little understanding of the 
underlying physical process. The investigation reported herein provides a more detailed 
insight into the problem to the extent of identifying the contributions of the various orders 
of scattering to the transfer of radiation through a slab. 

Apart from this, the investigation has certain practical connotations. In some of 
the more complex problems of radiative transfer involving (a) striated atmospheres in 
which scattering properties vary with depth, (b) scattering particles comparable in size 
with wavelength, and (c) polarization, methods of successive approximation may have to 
be utilized. In several of such methods (refs. 1 to 5), the successive iterations corre- 
spond to the contributions of successive orders of scattering. In applying these methods, 
knowledge of the role of multiple scattering would clearly aid in making a preliminary 
assessment of the rapidity of convergence of the procedure and hence the time and effort 
that would be needed to achieve a solution. Even in the case of simpler problems. 



sometimes it suffices to consider only the lowest orders of scattering (primary, sec- 
ondary, or tertiary). For these situations one can arrive at an answer of acceptable 
accuracy with a considerable saving of labor. Here also a knowledge of the role of 
multiple scattering would aid in deciding where the use of such approximate methods is 
warranted. 

In the case of the transfer of radiation through a homogeneous plane parallel atmo- 
sphere, even in the simplest of cases involving isotropic scattering, an analytical 
approach cannot be carried beyond secondary scattering (ref. 6), and in order to deter- 
mine the contributions of the higher modes of scattering, numerical procedures are 
required. Such approaches do not lend themselves to a systematic parametric analysis 
to determine the dependence of the different orders of scattering on slab thickness, albedo 
for single scattering, and asymmetry factor. In view of this fact, consideration has been 
limited to the one -dimensional radiative transfer problem. For this case exact solutions 
exist and the problem reduces to that of expressing the exact solutions as power series 
in the albedo for single scattering. The formulas so derived permit the determination of 
the contribution of any order of scattering for any slab thickness, for arbitrary albedo, 
and for arbitrary asymmetry parameter. 

The nature of the radiative transfer problem is such that the results of this one- 
dimensional analysis might be expected to have more than simply qualitative relevance. 

In this regard it is to be borne in mind that the two -beam approximation (refs. 7 and 8), 
which is tantamount to a one-dimensional approach, has proven to be of considerable 
value in providing results of acceptable accuracy to many three-dimensional problems of 
radiative transfer. 

The one-dimensional transfer problem has special significance in the case of high- 
energy corpuscular radiation where secondary production is in predominantly near forward 
angles and transfer calculations are made as if all secondaries proceed in the forward 
direction (straightahead approximation). (See ref. 9.) The multiple-scattering series 
performs a special function in this case since it reduces the usual integro -differential 
transfer equation into a set of coupled differential equations. (See refs. 10 and 11.) 
Particular advantage of the multiple -scattering series has been noted in application to 
charged-particle transfer calculations where, because of ionization energy loss, the series 
has nearly converged after only the double scattering term. (See ref. 11.) This rapid 
convergence is in spite of the high ’"albedo” (multiplicities near 3) in these scattering 
events. 

The author wishes to acknowledge his indebtedness to J. W. Wilson of the Langley 
Research Center for drawing his attention to the relevance of this research to the problem 
of high-energy radiation transfer. 
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SYMBOLS 


= f\t) + F^t) (proportional to photon density) 

F*(t) + F*(t) 

= (fraction of photons at depth t of nth order) 

F*(t) + f\t) 

= F*(t) - F^t) (net downward flux) 
upward flux of diffuse radiation 
downward flux of diffuse radiation 

contribution of nth order photons to reflected flux divided by incident flux 


D(t) 

Dn(T) 

F(t) 

F»(t) 

f‘(t) 

Fi(0) 

F^(t),Fj^( 7) contributions of nth order photons to upward and downward flux divided 
by incident flux, respectively 


'n(-o) 


F^(0) 


contribution of nth order photons to transmitted flux divided by incident 
flux 

contribution of nth order photons to reflected flux divided by total reflected 
flux 


FjI(to) contribution of nth order photons to transmitted flux divided by total 
transmitted flux 


g 






Il(r) 

l2(r) 


incident flux 

asymmetry parameter 

specific intensity of diffuse radiation 
»27T 


= r i(T,/i,0) d0 

specific intensity at depth r in direction of incident radiation 
specific intensity at depth t in direction opposite to incident radiation 


3 



k 




1/2 


p{li,(p;ii' ,(p') phase function 








^0 


k - 1 + CO 
k + 1 - CO 


0 

0 


6 Dirac delta function 

9 angle to upward vertical 

II = cos 6 


Mo 

T 

^0 

0 

00 

COq 


cosine of angle of incident beam with downward vertical 
normal optical depth 

normal optical depth of homogeneous plane parallel atmosphere 
azimuthal angle 

azimuthal angle of incident beam 

albedo for single scattering (elastic scattering probability) 
Primes denote variables of integration. 


ANALYSIS 


In considering the contribution of various orders of scattering to the radiative 
transfer through a homogeneous slab, the following recursive approach is invariably 
adopted. From the intensity distribution of the incident radiation, one calculates the 
intensity distribution of the first-order scattered photons. With this as an input, one 
then determines the intensity distribution of second-order scattered photons and so on. 
When one endeavors to apply this procedure in practice, problems arise. Thus, in the 
transfer through a plane parallel atmosphere, even in the simple case of isotropic 
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scattering, successive integrals become analsrtically unmanageable beyond the second 
order. (See ref. 6.) Even in the one-dimensional problem, although the successive 
integrations involve only elementary functions, the expressions become excessively com- 
plicated after a few terms. In view of this complication, an alternative approach has been 
adopted in the present paper. For a particular (albeit restricted) class of phase functions, 
the radiative transfer equations can be solved exactly. The problem then reduces to the 
mathematical one of developing the exact solutions as power series in the particle albedo. 
The successive terms of such a series development give, of covurse, the contributions of 
successive orders of scattering. The more detailed development of this line of attack 
is pursued herein. 

The following equation for the transfer of diffuse radiation through a homogeneous 
slab illuminated from above by a parallel beam making an angle cos“l with the 
downward vertical is derived in reference 12: 

1 2tt 

^dl(j^x^^ - dfX’ r p(/i,<|>;|Li>') I(t,M', 0’) d<^)' 

dr 4n JQ 






( 1 ) 


In this equation cos“l /j. is measured from the upward vertical and the phase 
function is normalized in the sense that 

^ 2 tt pi 

J ^ diUdc^) = 47ro)Q 

Integrating equation (1) over <j) yields 


ttF 


0 (0), V ~V>^0 

-p^ 


where 


= r 1{t, II, < p) d0 

do 


( 2 ) 
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and 


Bear in mind that the phase function depends only on the cosine of the angle of scattering 
<p - 

One now fvirther restricts the phase function to the form 

= WqQi + g)6(ju’ - p) + (1 - g)6(M' + /i)) (3) 

where cOq is the albedo for single scattering.^ 

By substituting this expression for p^®^(jLt,p') into equation (2), one obtains 

M - -^|(1 + g) + (1 - g) 

("^0 ~ ■ S)6(-Mo + e (4) 

It is to be noted that the forcing function (last term of eq. (4)) Involves Dirac delta func- 
tions, and this same characteristic must perforce be imposed on the solutions. Denoting 
the intensity in the direction of the incident radiation by and in the opposite direc- 

tion by I2 (t) yields 

= I^(t) 6(|u - Pq) (p > 0) 

= l2(r) 6 (m - Mo) (M > 0) 

By using equations (4) and (5), the equations giving the downward and upward inten- 
sities are, respectively: 

^In the "one-speed approximation" in neutron transport theory, coq would be the 
average multiplicity. 
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dIi(T) 

- Mq) = - Mq) ■ ■*■ ■ ^o) ■*■ ^^ ■ S)l2('^) 6(M - Mq)^ 


77FoWq 


(1 + g)5 (ii - Mq) + (1 - g)5 (-/^Q - m) 




( 6 ) 


dloCr) ojq 

M — ^ - /^o) " '^o) " Y’ " f^o) + ■ g)Il('^) ®(M - Mq) 


TTFoC^or 


(1 + g)6 (-iIq - + (1 - g)6 ^-iiq + fxj 


-t/Mq 


(7) 


Integrating over their respective hemispheres yields 


Mo 


dl.(T) 


tiTi ^nT *1 “VMn 

^ = -Ii(t) + ^[(1 + g)Ii(r) + (1 - g)l2(T)J + — ^ (1 + g)e ^ ° 


( 8 ) 


and 


dlgCr) Wf,r TTF^Wn 

Mo = l2^^) - (1 


g)e 


-t/Mo 


(9) 


By introducing downward and upward fluxes as defined by 
MoIjCt) = 7tF*(t) 


and 


Mol2(’’) = ^F^(t) 


into equations (8) and (9), one arrives finally at the pair of first-order linear differential 
equations 


dF‘_ + 

'"o-d7^ ^ 


(1 +g)F‘ + (l - g)F*j+- Q^Q‘^^ (1 +g)e 


( 10 ) 
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and 


Mq ^ = F* - ^[(1 + g)F* + (1 - g)F j (1 _g)e 


-t/Mo 


In addition, the following boundary conditions are to be imposed: 


( 11 ) 


o 

ll 

(12) 

f'(To) = 0 

(13) 


Since the postulated phase function as given by equation (3) requires that all photons 
be scattered in either the forward or the backward direction, not unexpectedly equa- 
tions (10) to (13) are precisely those governing one -dimensional radiative transfer. Exact 
solutions for these equations have been derived by Sobolev (ref. 13). For the convenience 
of the reader, these solutions are rederived in appendix A. 

It is of interest to note that the phase function as given by equation (3) with g « 1 
approximates Mie scattering from large smooth spheres and with g = -1 scattering 
from opaque particles with rough faces oriented at random (ref. 14). 

Throughout the subsequent analysis, Pq is set equal to unity (that is, normal 
incidence). However, by simply replacing t by t//j.q in any of the results derived 
below, one obtains answers corresponding to the case in which the slab is obliquely 
illuminated at an angle cos~l Pq to the downward vertical. 

In appendix A the following expressions are obtained for the overall intensities of 
the downward and upward directed beams at a depth t within a slab of optical thick- 
ness Tq, particle albedo Wq, and asymmetry factor g: 





e 


-T 


(14) 


and 



^k(ro-r) -k(r„-T) 


kTfi 


2 

- V® 


(15) 
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where 


and 


( 16 ) 


k - 1 + 0 )q 
® k + 1 - o)q 


(17) 


Contribution of Successive Orders of Scattering 
to Reflection From a Semi-Infinite Slab 

By setting Tq = in equations (14) and (15), the following expressions for the 
intensities as a function of optical depth are obtained: 


and 


F*(t) 


= e 


-kr 


(18) 


F‘(t) _ ,, ^-kT 


= roe 


(19) 


The reflected flux, divided by the fltix of the incident beam, is thus 

f\o) _ 

^0 

By using the series development of rg as given in equation (B20), one obtains 

n n+1 




(2n - 2j)! 


(n - j)'.(n - 2j + 1)1 


4g 


(1 + g)' 


( 20 ) 


Hence, the contribution of nth order scattering to the reflected beam (normalized by 
dividing by the incident flux) is given by 


9 



( 21 ) 




n n+1 
2’ 2 

V (2n - 2i)! 

4g "j 

4 J 

' li-gj 

j=0 j’(n - j)’.(n - 2j + 1)! 

^ (1 + g)?J 


By setting n = 1, 2, and 3, one obtains the following contributions of the first three orders 
of scatterii^: 

First order: 

-^(1 “ g) (22) 


Second order: 



(1 + g)(i - g) 


(23) 


Third order: 




(l+g)2(l-g) + i(l 



(24) 


These expressions not only define the contribution of successive orders of scattering but, 
in addition, reveal something of the nature of the scattering processes that are involved. 
Thus, the presence of the factor 1 - g in expression (22) is a mathematical manifesta- 
tion of the rather obvious fact that the first-order photons in the reflected beam have all 
undergone a single rearward scattering. In the case of the second-order photons, the 
presence of the factor (1 + g)(l - g) in expression (23) indicates that these have all 
undergone one forward and one rearward scattering (although not necessarily in that 
order). The third -order photons in the reflected beam are of two kinds correspondii^ to 
the two terms of expression (24). The first kind embraces those which have undergone 
two forward scatterii^s and one rearward scattering and the second kind, those which 
have been subjected to three rearward scatterings. Note that in the emerging beam the 
former preponderates over the latter by a factor of four. 


Contribution of Different Orders of Scattering to Photon 

Population at any Depth Within a Semi -Infinite Slab 

The density of photons (D) at a depth t is proportional to F*(t) + F*(t). 
from equations (18) and (19), 


Thus, 
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D = Constant 




However, from equations (B14) and (B20) 

f A„V') 


n=l 


and 


•o = Z ®n"o" 


where 


An = 


(1 +g) 

2 


n=l 


vn n r 


n-1 n 
2 ’2 


(25) 


(26) 


(27) 


n y ^ y (-)" . y 


m( 


. . Cm - 2(n - 3 - 1)3 

r 7 

iUn - 2j)'. 

11 + g)^ 


(28) 


n n+1 
2 ’ 2 


B =lLt3.t(l±A\ y (2n - 2j)! r 4g _~ ) 

" \ 4 / \1 - gj j'.(n - j)’.(n - 2j + 1)1 L (1 + g)^ 


(29) 


Substituting equations (26) and (27) into equation (25) yields 


D = Constant 




.-T 


n=l 


n=l 


n-1 


= Constant ^ e~^ + ®n + ^ A.B^.j COq" 

n=l V 3=1 
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Dividing by equation (25) yields 


1 - 

n=l 


( 30 ) 


where Dj^ 
given by 


defines the fraction of photons at depth 


D. 


n-1 

■^n + ®n ^ ^j^n-j 



(l.ro)e(l-k)--l 


T which are of nth order and is 


(31) 


Contribution of Various Orders of Scattering to Reflection From 


and Transmission Through a Finite Slab 

The reflected flux (normalized with respect to the incident flux) is obtained on 
setting T = 0 in equation (15) 



-2kr 

1 - e 0 
1 - rp2 e C 


(32) 


The transmitted flux (normalized as above) is obtained on setting t = Tq in 
equation (14) 


F 



F 


0 



(33) 


Consider the case Wq = 1 (that is, conservative or pure elastic scattering). From 
equations (16) and (17) k = 0 and Tq = 1 and the expressions (32) and (33) for the 
reflected and transmitted intensities, assume the indeterminate form 0/0. 

In this case let Wq = 1 - e where e is supposed to be sufficiently small to justify 
the neglect of all terms but those involving the lowest order in e. Thus, 


k = \/l - ge 


1/2 


r. = l - -^el/2 

- g 


(34) 
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The introduction of these expressions into equation (32) yields 



- g)'To 

2 + (1 - g)TQ 


(35) 


The introduction of equations (34) into equation (33) yields 



- ^ - , 1/2 

\/l - g 

— ^ e^/2 + 2^jl - g e^/2 Tq 

- g 


2 

2 + (1 - g)TQ 


e 


-T 


0 


(36) 


By the substitution of equations (26) and (27) into equations (32) and (33), one can, by 
straightforward procedures similar to those used in the preceding subsection, express 


ZlM ai.d — 


0 


'0 


as a power series in Wq. Thus, 


= ! <(«) V 


0 n=l 


(37) 
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( 38 ) 



By dividing equation (37) by equation (32) ^ 1^ or by equation (35) (Uq = 1), one 

obtains 

OO 

1 = ^ F>) 

n=l 


A n \ / 

where F '(0) = — - — and thus defines the fractional contribution of the nth order photons 
F‘(0) 

to the reflected beam. 

Similarly, by dividing equation (38) by equation (33) ^for u)q ^ 1^ or by equa- 
tion (36) ^for Wq = 1^, one obtains 





n 


where 


\ 

n 



llM 


and thus defines the fractional contribution of nth order photons to the transmitted beam. 


NUMERICAL RESULTS 


The formulas of the preceding section have been used to make calculations bearing 
on the role played by the various orders of scattering on reflection from semi-infinite 
slabs and reflection from, and transmission through, slabs of finite thickness. In this 
section the results of these calculations are presented and discussed. The effects of 
variation in slab thickness, particle albedo, and asymmetry factor are considered. 

It is of interest to note that, in practice, both the albedo for single scattering Wq 
and the asymmetry factor g are subject to wide variation. Thus, the albedo for single 
scattering can be as high as 0.9999 in clouds. Insofar as the asymmetry parameter is 
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concerned, this value approaches unity for Mie scattering from large smooth spheres 
and minus unity for opaque particles with rough faces oriented at random. 

Overall Reflection From Semi-Infinite Slab 

In figure 1 the reflection coefficient is plotted against the asymmetry factor for 
various values of particle albedo coq. As the asymmetry factor increases (the 
scattering increasing in the forward direction), the impinging photons penetrate deeper 
and, thus, for particle albedos other than unity, the absorption increases and the reflec- 
tion decreases. When the scattering is entirely in the forward direction (g = 1), the 
reflected flux must, of course, fall to zero. At a particle albedo of unity, there is no 
absorption and, since the steady case is being treated and there is no accumulation of 
photons within the slab interior, there must be as many photons emerging as are entering. 
In other words, the reflection coefficient is unity independent of the value assigned to the 
asymmetry factor. In view of this effect, one would expect the reflection to become 
insensitive to asymmetry factor at high particle albedos. However, it is interesting to 
note how large an albedo is needed to produce a reasonable measure of insensitivity to 
the asymmetry factor. Thus, even with a particle albedo as high as 0.99, the reflected 
flux at g = 0.8 is only about 75 percent of that at g = -0.8. 

Overall Reflection From and Transmission Through Finite Slabs 

For the case of conservative (elastic) scattering ^coq = 1^, as many photons leave 
the slab as enter it, since the time-independent problem is being considered. Thus, 

Fq = F^e'^O + f‘(0) + F*(to) 

The first term appearing on the right-hand side corresponds to those photons which pass 
through the slab without undergoing a single scattering. The remaining two terms corre- 
spond to those photons which are diffusively reflected from and transmitted through the 
slab, respectively. In the case of thick slabs, the first term becomes negligibly small 
and the preceding relation then becomes 

1 F*(0) , 

^0 ^0 

Figure 2 shows dependence of reflection from and transmission through a slab of optical 
thickness 0.1 on asymmetry factor for values of particle albedo of 0.2 (0.2) 1.0. Similar 
plots for slabs of optical thickness 2 and 20 are given in figures 3 and 4, respectively. 
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In figure 4(a), note how sensitive the reflection coefficient from a thick slab 
becomes to particle albedo as conservative scattering = Ij is approached. Thus, 
for isotropic scattering (g = 0), an increase of particle albedo from 0.9 to 1.0 results in 
almost a doubling of the intensity of the reflected beam. Transmission through a slab 
of optical thickness 20 only becomes significant at particle albedos in excess of 0.9. 

Note how sensitive transmission is to particle albedo at these high values. Thus, again 
considering the case of isotropic scattering, an increase of particle albedo from 0.99 
to 1.0 results in a doubling of the intensity of the transmitted beam. Note also the sensi- 
tivity of transmission to asymmetry factor in the neighborhood of g = 1. In the case of 
conservative scattering, an increase in the asymmetry factor from 0.8 to 1.0 results in a 
threefold increase in transmission. 

Contribution of the First 10 Orders of Scattering to 
Reflection From a Semi-Infinite Slab 

The contribution of the first 10 orders of scattering to reflection from a semi- 
infinite slab is plotted against particle albedo over the range 0.8 to 1.0 for various values 
of the asymmetry factor in figure 5. Orders of scattering in excess of 10 only play a 
significant role at the higher albedos. Thus, at a particle albedo of 0.8, the first 10 orders 
of scattering contribute more than 90 percent to the reflected beam for all values of the 
asymmetry factor. However, at an albedo of unity and predominantly forward scattering 
(g = 0.9), photons which have been scattered more than 10 times contribute 80 percent 
to the reflected beam. 

Contribution of the First 10 Orders of Scattering to Reflection From 
and Transmission Through Slabs of Finite Thickness 

The contribution of the first 10 orders of scattering to the reflected and trans- 
mitted beams is given as a function of slab thickness for various values of particle 
albedo in figure 6 for the case of predominantly forward scattering (g = 0.8). One notes 
from figure 6(a) that insofar as the reflected beam is concerned, the modes of scattering 
in excess of 10 only play a significant role at particle albedos in excess of 0.8 and slab 
thicknesses in excess of 2. Indeed, for sufficiently thin slabs, only first-order scat- 
tering will contribute. For a particle albedo of 0.9 as the slab thickness increases 
beyond 2, orders of scattering in excess of 10 make an increasing contribution. The con- 
tribution tends to a limiting value as a slab thickness of about 14 is approached. At this 
thickness, the reflected beam is, to all intents and purposes, the same as that for a semi- 
infinite slab. At a particle albedo of unity, and a slab thickness of 20, photons that have 
undergone more than 10 scatterings contribute 50 percent to the reflected beam. Further 
increase of slab thickness will lead to further enhancement of the contribution of these 
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higher modes of scattering. By referring to figure 5, it is to be noted that the limiting 
contribution in this case is about 66 percent. 

By turning to a consideration of transmission (fig. 6(b)), one notes that the orders 
of scattering in excess of 10 only assume a significant role for slab thicknesses in excess 
of 2. As the slab thickness increases, there will be, of course, an exponential fall off in 
the transmitted flux. However, of those photons that do emerge, the higher orders of 
scattering will make an increasingly important contribution. Thus, the curves, even for 
small particle albedos, will tend monotonically to zero. This is simply a manifestation 
of the fact that although for small particle albedos the probability of a photon surviving 
more than 10 scatterings is low, for sufficiently thick slabs the probability of photons 
getting through at all with 10 or fewer scatterings is even smaller. 

Similar plots are presented for the case of predominantly rearward scattering 
(g = -0.8) in figure 7. Comparison of figures 6(a) and 7(a) and 6(b) and 7(b) indicated, 
by and large, the more important role of the higher orders of multiple scattering associ- 
ated with strongly forward scattering. In figure 8 the contributions from orders of 
scattering in excess of 1 to reflection from and transmission through a slab of optical 
thickness 0.1 are plotted against particle albedo for the cases of predominantly forward 
(g = 0.9) and rearward (g = -0.9) scattering. Even for this thin slab for g = -0.9 and a 
particle albedo in excess of 0.8, orders of scattering of 2 or higher contribute over 
40 percent to the transmitted beam. Similar plots are presented in figure 9 for a slab 
of optical thickness 20. In this case the contributions from orders of scattering in excess 
of 10 are displayed. It is interesting to note that at a particle albedo of unity and with 
dominantly rearward scattering, orders of scattering in excess of 10 contribute more 
than 99.99 percent to the transmitted beam. 

Dominance of Various Orders of Scattering With 
Depth in a Semi-Infinite Slab 

In the photon population at any depth, all orders of scattering are presented. As 
the depth increases, the principal contribution will come from increasingly high orders 
of scattering. This is exhibited in figure 10 for the case of isotropic scattering in which 
the domains of dominance of the various orders of scattering are charted within a frame- 
work of particle albedo and optical depth. Thus, for particle albedos less than 0.2, first- 
order scattered photons preponderate down to a depth of 6.6. Between 6.6 and 9.4, 
second-order scattered photons provide the dominant contribution and so on. As the 
particle albedo increases, the transition from dominance by one order of scattering to 
the next as the depth increases becomes increasingly rapid. Thus, at a particle albedo 
of unity, first-order scattering makes the major contribution down to a depth of 2.7; 
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however, by the time a depth of 7.5 is reached, the greatest contribution comes from 
tenth order of scattering. 

Dominance of Various Orders of Scattering Within 
Beams Transmitted Through Finite Slabs 

As the slab thickness increases, the principal contribution to the transmitted beam 
will come from increasingly high orders of scattering. This is illustrated in figure 11 
for cases of predominantly forward and rearward scattering, respectively. In each 
instance the domains of dominance by various orders of scattering are displayed within 
a framework of particle albedo and slab thickness. In figure 11(b) one notes that although 
first-order scattering plays a major role for small particle albedos and small slab thick- 
nesses, thereafter dominance is invariably by an even order of scattering. This is a 
manifestation of the fact that for predominantly rearward scatterir^, the odd orders of 
scattering are suppressed in the transmitted beam. Further manifestation of this effect 
is found in the plots described in the next subsection. 

Contribution of the Various Orders of Scattering to Reflection From 
and Transmission Through Slabs of Finite Thickness 

In figure 12, plots are given that show the individual contributions of the first three 
orders of scattering to reflection from and transmission through a slab of optical thick- 
ness 0.1 for albedos of 0.2 and 1.0 with predominantly forward scattering. For this thin 
slab the reflected and transmitted fluxes consist almost entirely of first-order photons. 
From figure 13 which corresponds to predominantly rearward scattering, one notes that 
at high particle albedos second-order photons contribute as much as 30 percent to the 
transmitted beam. 

Figure 14 displays the contributions of the first 10 orders of scattering to the 
reflection from and transmission through a slab of optical thickness 20 for various par- 
ticle albedos for the case of predominantly forward scattering (g = 0.8). Insofar as 
reflection is concerned, the contributions of successive orders are monotonically 
decreasing. However, as the particle albedo increases, the spectrum becomes pro- 
nouncedly flatter and subscribes once again to the increasingly important role played by 
all higher modes of scattering at the higher particle albedos. The curves pertaining to 
transmission, on the other hand, are peaked. Thus, at a particle albedo of 0.2, third- 
order scattering contributes most. At a particle albedo of 0.4, it is the seventh order of 
scattering which preponderates. At particle albedos of 0.6 and higher, the principal con- 
tribution comes from orders of scattering in excess of 10. 
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In figure 15 similar plots are presented for the case of dominantly rearward scat- 
tering. In this case it is to be noted that the curves appear to have a jagged or saw- 
toothed structure. This is a further manifestation of the fact that with predominantly 
rearward scattering, the even orders of scattering are suppressed in the reflected beam 
and the odd orders in the transmitted beam. This is strikingly illustrated in figure 16 
in which plots are given for a slab of optical thickness 2 and an asymmetry factor of 
g = -0.98. 


Orders of Scattering Required for Specified Accuracy 

Figure 17 provides plots that define the order of scattering needed to give an accu- 
racy of 1 percent in the calculation of reflected and transmitted fluxes for a slab of opti- 
cal thickness 0.1. For reflection with predominantly forward scattering and transmission 
with predominantly rearward scattering, second-order terms must be included at all but 
the smallest albedos if the desired accuracy is to be achieved. 

Similar plots, defining order of scattering needed for 10-percent accuracy, are 
given in figures 18 and 19 for slabs of optical thickness 2 and 20, respectively. Insofar 
as reflected flux from slab of optical thickness 2 is concerned, the case of predominantly 
forward scattering is the severest one since higher orders of scattering are needed to 
achieve the desired precision. Thus, at a particle albedo as low as 0.2, second order must 
be included and at high particle albedos, up to fifth order must be considered. Insofar as 
the transmitted flux is concerned (fig. 18(b)), the roles are reversed and it is the case of 
predominantly rearward scattering which requires the inclusion of the higher orders of 
scattering. Very high orders of scattering must be embraced in the case of the thick 
slab ^Tq = 20^ if the desired precision is to be achieved. For g = 0.9, orders of scat- 
tering in excess of 10 are needed for a particle albedo exceeding 0.8 in the computation 
of reflection and for a particle albedo exceeding 0.33 in the computation of transmission. 

CONCLUDING REMARKS 

In the present paper expressions for the contributions of various orders of scat- 
tering to the overall reflection from and transmission through slabs of arbitrary thick- 
ness, particle albedo, and asymmetry parameter have been given for the one -dimensional 
case and numerical examples have been calculated. 

The results presented herein have several applications. In the first place, they 
will aid in identifying those physical situations - combinations of thickness, albedo, and 
scattering phase function — in which answers of adequate accuracy can be achieved by 
considering only first or second orders of scattering. If this is known, solutions of 
problems in radiative transfer can be obtained with relative ease. 
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The second area of application relates to those situations where the complexity of 
the problem necessitates the use of numerical procedures and where the procedure 
adopted is such that successive iterations define the contributions of the successive 
orders of scattering. The results contained herein will permit a preliminary assess- 
ment to be made of the rapidity of convergence of such series developments and, hence, 
the effort that will be needed in achieving a solution of the desired precision. 

A third area in which methods similar to those discussed herein may find applica- 
tion is to high-energy radiation transfer in which the high-energy secondaries are pro- 
duced predominantly forward. However, low-energy secondaries are produced and sub- 
sequently scatter nearly isotropically. The approximation of low-energy propagation by 
means of a one -dimensional isotropic scattering law (asymmetry factor zero) would 
surely produce more realistic results than a strict straightahead approximation 
(asymmetry factor unity). 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
October 2, 1975 
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APPENDIX A 


EXACT SOLUTION OF RADIATIVE TRANSFER EQUATION 
IN ONE DIMENSION 

It has been shown in the text that by postulating a phase function of the form given 


by equation (3), that is, 

p^®\p,]u’) = - m) + (1 - g)5(M' + m3 

the equations of radiative transfer in the case of normal illumination given 

by the following equations (eqs. (10) to (13)): 

= .F ' + * g)F ' + (1 - g)F^+ 1^(1 + g)e-^ (A2) 

dr 2 L J ^ 

111! = F ♦ - (1 + g)F* + (1 - g)F * ^-^(1 - g)e”^ (A3) 

F*(0) = 0 (A4) 

F<(tq) = 0 (A5) 

The addition and subtraction of equations (A2) and (A3) gives 

|-(f» + f‘) = -^1 - gcoQ^ (f» - F<) + FQgujQe"^ (A6) 

f') = -(1 - (f' + f')+ Foa,„e-' (A7) 


Denote F* + F* (proportional to energy or photon density) by D and F* - F* 

(net downward flux) by F and introduce these quantities into equations (A6) and (A7). 
The results are 
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^+(l-go.o)F = Fog»oe-" 


(AS) 




(A9) 


These two simultaneous linear equations with constant coefficients can be solved by 
using standard procedures. A particular solution is, obviously, 




D^-F^e-r 


F = -Foe-^ ^ 


(AlO) 


By turning to the homogeneous equations, one finds 



Try 


D = de^'^ 

F = fe^'^ 

Hence 

kd + - gWQ^f = 0 

^1 - ojQ^d + kf = 0 
These relations require that 

k2 = (l - cOq) (l - g^o) 

= 1 - (1 + g)WQ + g(^g2 


(All) 

(A12) 


(A13) 
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and the general homogeneous solution takes the form 


D = d,e"^^ + doe^"^ 


F = f + f2©^^ 


In order to establish the interrelationships between these four constants, it suffices to 
substitute them into equation (A12) 


(l - coqJ (dje-kT + dgekr) . k(fj 


e-kr _ f„ekT\ = o 


'Renee 


, ‘'"Oh 


1 - Wq 

f, d« 

2 k 2 


By substituting these equations into equations (A14) and adding the particular inte- 
gral, one finds 


D= F*+ F* = die"*^"^ + doe*^T^ - F^e 


kr ^ H okr 


F = F * - F* = 


1 - 0 ) 


2(d,e->'^ - - Fq®' 


The addition and subtraction of equations (A15) and (A16) gives 


, /k + 1 - cOfA , /k - 1 + (jOf\ , 


- FQe”’’’ 
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F»(r) = dj 


- 1 + cOf\ . /k + 1 - (jj(\ - 


(A18) 


The two remaining constants are determined by enforcing the bovindary condi- 
tions (A4) and (A5), 


'k + 1 - o)j 

2k 


+ dr 


'k - 1 + coj 
k 2k 


-Fo = 0 


(A19) 




2k 


(A20) 


From equation (A20) 


-2kT- 

d2 = -djrQe 0 


(A21) 


where 


^0 = 


k-l+Wj- 
k + 1 - Wr 


(A22) 


Substituting equations (A21) and (A22) into equation (A19) yields 


'k + 1 - coj- 

I 2k 


1 - r 

To e 


Hence, 


rk + 1 - 0 )(- 
I 2k 


-2kr 

2 ■^’'0 ° 
1 - rQ^^e ^ 


By inserting these constants into equations (A17) and (A18), one finally obtains for the 
downward and upward fluxes 
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and 


f‘(t) _ 


- Tq e 


kTf, „ 

e ^ - ^0^^ 


-krn 


- e 


-T 


fHt) _ 
^0 


° e ''"0 . ,^ 2 e -‘‘"0 


(A23) 


(A24) 
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EXPANSION OF e ^ ' AND AS POWER SERIES IN Wq 

The functions and are defined in equations (A13) and (A22) as 

follows: 


and 


r 91^/2 

k ( wo ) = J 


(Bl) 


k - 1 + 


'o(“o) = f7TT 




(B2) 


- k (" o ) 


as a Maclaurin Series in Wf 


Expansion of e its it iVlcU^ldUX 111 iJCXtCS ill UL/q 

In the following derivation, use is made of the following expression for the nth dif- 
ferential coefficient of a function of a function. 


If y = f(u) and u = , then 


n'- dWn ” r’. ^ 


0 r=l 


where is the coefficient of h^ in 


nr 


0(o)Q -hh) - 


In the present case 

f(u) = f^^^(u) = e^ 


(B3) 


(B4) 


where 


u = - ^1 - Awq - Bo)Q^j 
A = 1 -I- g 


1/2 


B = -g 


(B5) 
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Determination of The coefficient ^Kj, appearing in equation (B3) may be 

determined as follows: 


0(^0 +h) - 0 (wq) 


= ^^(i-ao>o-bV) 


1/2 


r 2]1/A 

- 1 - A(cJq + h) - B(coQ + h) J J 


= T 



1 — — 2)^^^ 

(1 - 1 

1 - (l - Ah - Bh^j 


(B6) 


where 


A + 2o;qB 

1 - Awq - Bcoq^ 


A = 


B= B 

1 - Awq - Bojq^ 

The binomial theorem states that 


(B7) 


(1 +x)^= 1 + y r(r - 1) . . . (r - p + l l^p 

^ p'. 

p=l 


and when r is an integer, the series is a terminating one and is given by 


( 1 +xf = f r! ^p 
L/_ pt (r - p)! 


p=0 


(B8) 


(B9) 


By using equation (B9) in conjunction with equation (B6), one finds 
0(o.o +h) - 0(0^0^ = tt(i - Acoq - Bcoo 2^ ^ . Ah - Bh^) 


m=0 


(BIO) 
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Applying the binomial theorem for nonintegral exponent (eq. (B8)), one finds 

, -.m/2 V ('f) (^ ■ f) • • ■ ^ - 2 \^ 

(l - Ah - Bh2) =1+2 ■■ ^ 

p=l 


Making use of equation (B9) results in 


^^m/2 ^ (■t)(^ ‘ y) ■ ■ ■ ^ ” y)/- 'iP V p' 

(l - Ah - Bh2) = 1 . ^ A ^(Ah) 2 

p=l q=0 


n-1 n 


, V M-?)l--f) (■--'- 8.-.,., 

' A A j’.(n-2j)! 


n=l j=0 


(Bll) 


1 Y\ 

where the upper bound is ~ or - whichever is an integer. 

/ _ _ o\m/2 „ 

Thus, in the expansion of (1 - Ah - Bh'^j , the coefficient of h“(n ^ 1) is 

given by 

n-1 n 
2 


j’. (n - 2j): 


j=0 


(B12) 


By bearing in mind that is the coefficient of h” in <p^Q + h^ - 

obtains from equation (BIO) and expression (B12) 




, one 
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n-1 n 


! y (-)“ 



Lj m! (r - 

m)' l< j'. (n 

- 2j)'. 


m=0 

]=0 



n-1 n 




^ 2 ’2 

(-)™ r'. V m(m - 2) . 

. . . Qn - 2(n - ; 

j - li) / 2 bV 

0 y Z/ 

m=0 

m! (r - m)l Lj 

j=o 

j'. (n - 2j)'. 

U 2 ; 


(B13) 


Maclaurln expansion of 
into equation (B3) yields 



T 

The substitution of equations (B4) and (B13) 


/2 H 

^ d"y ^ y y r: \ ^ - 2) . . . Qn - 2(n - j - 1)) / 2 bV 

n- dwr" V'2/ r\ mil 


r=l 


£i i 

^ V 

- m)'. ^ 


— „ ml (r - m)'. — „ 

m=0 ' j=0 


j'.(n - 2j)'. 




When cOq = 0, then k = 1, u = -r; A = A = 1 + g and B = B = -g 


n-1 n 



-TTd +g)l 

” V ^ V (-)™ m(n'> - 2) . . . (m - 2(n - j - 1)} 

r 2g 1 

L 2 J 

kx U J’(n-2i)'. 

ji + g)?J 


Thus, 




(B14) 


where 


n-1 n 


r (1 + g) 

n n r 2 ’2 , 

V ^ V <-S^ \ m(m - 2) . . . 1 

m - 2(n - j - l'^ 

r 2g 1 

L 2 J 

4 Zo ='<" 

- 2j)! 

.U + g)j 


(B15) 
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Expansion of Tq^coq^ as a Power Series in u)q 
O ne may rewrite equation (B2) in the form 

_ k - (1 - (Op) 

"■O " k + (1 - c^o) 

k2 + (1 . «,„)2 - 2k(l - «,) 
k2-(l-c.o)2" 

However, from equation (A13) 

k2 = (1 - o^o) (1 - gu>o) 

By substituting equation (B17) into equation (B16), one obtains 


(B16) 


(B17) 


^0 = 


l-fifEcp-k 


1 - g 


Y~i^o 


(B18) 


Setting m = 1 in equation (Bll) one obtains 


k = 1 1 - Acoq - Bcoq 


n-1 n 


1/2 


= 1 + J y ' t: ‘ ' 'il A"-2iBi 

A ® A j- (n - 2j)! 


n=l j=0 


n-1 n 


n=l 



n ^ ^ 

1 V (2n - 2j - 2)! 1 

+ g\ 

\2j 

Af. j’-(n - 2j)!(n - j - 1)! 

3=0 

1 2 ) 


\n-2j 


(B19) 
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The substitution of equation (B19) into equation (B18) yields 


n-1 n 
n 2 ’2 


^0 = 


/ \ ) ^ \ ^ ^ 


(2n - 2j - 2)1 


(n - 2j)t (n - j - 1): \ 2 






n-1 n 

n 2 ’2 


y y (2n-2j -2)! fl -»■ _g)3 

L\2) L jI(n-2j)!(n-j-l)lV 2 ) ^ 

n=2 ]=0 


The adoption of s = n - 1 as the running variable gives 


o=!(h 

S=1 ' 


S S+1 


^ »„Y 1 

fl +g'^ 

1 V (2s - 2j)! 

r 1 

4 0 } ' 


' j'(s - j)'.(s - 2j + 1)! 

1=0 

^(1 +g)?. 


(B20) 


Case of Isotropic Scattering (g = 0) 

In this case the only nonvanishing term in the second summation corresponds to 
j = 0 and 


r 


0 


y ^ (2s)I — 

b 22ss.(s + i)t 


(B21) 
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Figure 3.- Reflection and transmission coefficients for a slab of optical thickness 2 as a function 

of asymmetry factor and particle albedo. 
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Figure 4.- Reflection and transmission coefficients for a slab of optical thickness 20 as a function 

of asymmetry factor and particle albedo. 
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Figure 6.- Contribution of first 10 orders of scattering to reflection from and transmission 

through finite slabs, (g = 0.8.) 
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Figure 7.- Contribution of first 10 orders of scattering to reflection from and transmission 

through finite slabs, (g = -0.8.) 
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Figure 8.- Contribution of photons that have undergone more than one scattering to reflection and 

transmission for a slab of optical thickness 0.1. 
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Figure 11.- Term in series expansion contributing most to transmitted beam. 
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Figure 12.- Contributions of successive orders of scattering to reflection and transmission for 

a slab of optical thickness 0.1. (g = 0.8.) 
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Figure 14.- Contribution of successive orders of scattering to reflection and transmission for 

a slab of optical thickness 20. (g = 0.8.) 
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Figure 15.- Contribution of successive orders of scattering to reflection and transmission for 

a slab of optical thicloiess 20. (g = -0.8.) 
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Figure 17.^ Order of scattering needed to give within 1-percent accuracy. 
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(b) Transmission: optical thickness 2. 

Figure 18.- Order of scattering needed to give within 10-percent accuracy. 
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Figure 19.- Order of scattering needed to give within 10-percent accuracy. 
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